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Abstract. The accuracy in position and flux density of 19 large radio continuum source catalogues has been 
determined using SPECFIND, a new tool recently made available through the CDS. The ~ 67 000 radio continuum 
spectra with three or more frequencies produced by SPECFIND were used to cross-correlate sources from different 
catalogues and to calculate offsets in right ascension and declination in the various catalogues with respect to the 
positions given in the NRAO VLA Sky Survey (NVSS) catalogue, which was adopted as a reference. The flux 
densities reported in the catalogues were compared to those predicted by the composite spectra, enabling us to 
assess the quality of the flux density calibration of the different catalogues. 
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1. Introduction 

A basic, fundamental source of information for the eval¬ 
uation and use of large radio continuum surveys is the 
position and flux density accuracy of the final source cat¬ 
alogues. The precisions in flux densities and positions can 
be determined through the examination of several sources 
of potential errors. We distinguish between direct errors 
connected to observations and errors due to the extrac¬ 
tion of sources from observed maps. The main sources of 
uncertainties are: 

— Position: telescope pointing errors, which can amount 
to < 10 % of the HPBW, as well as errors intro¬ 
duced by the algorithm for source extraction from 
the observed maps, i.e. usually the fitting of a two- 
dimensional Gaussian. 

— Flux density: errors in the absolute flux calibration 
(see, e.g., Ott et al. 1994), which are usually at the 
<10% level, as well as errors introduced by the 
source fitting algorithm - especially for extended, non- 
Gaussian sources, Gaussian fits can be insufficient or 
lead to the fitting of artificial multiple components; it 
is very difficult to obtain error estimates of the latter 
effects. 
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To determine the position accuracy of large radio 
source catalogues, these are in general cross-correlated 
with other published radio source catalogues, which have 
a higher resolution and which were generally obtained at 
different frequencies. Often used reference catalogues are 
the TeXas Survey of discrete radio sources (TXS, Douglas 
et al. 1996) and the Green Bank 4850 MHz survey (GB6, 
Gregory et al. 1996). To determine the flux density ac¬ 
curacy of catalogues, there are two possibilities: if there 
is another published catalogue at the same frequency, the 
flux density scales can be compared directly, and if there 
are (at least) two other suitable catalogues at different 
frequencies, the slope of the radio continuum spectra can 
be calculated based on these catalogues and the flux den¬ 
sity at the newly observed frequency can then be inter- or 
extrapolated. 

In a recent paper (Vollmer et al. 2004), we presented 
the new SPECFIND tool for the extraction of cross¬ 
identifications and radio continuum spectra from radio 
source catalogues contained in the VizieR database of the 
Centre de Donnees astronomiques de Strasbourg (CDS). 
The SPECFIND code can handle radio source catalogues 
made at different frequencies with different resolutions. 
For the radio spectrum a power law is assumed, which 
can have one break if more than two independent fre¬ 
quency points contribute to each of the two slopes. Since 
our cross-identification is based on the comparison of flux 
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densities at the same frequency and on the radio spec¬ 
trum at different frequencies, SPECFIND takes into ac¬ 
count whether a source is resolved by a given survey or 
not - if, e.g., a source is resolved only by our reference 
survey, the NRAO VLA Sky Survey (NVSS; Condon et 
al. 1998), and not by other available surveys, there will be 
no positive cross-identification, because the flux density of 
the resolved NVSS source cannot be fitted into the radio 
spectrum of the unresolved source as defined by the other 
catalogues. 

With this tool we identified ~ 67 000 radio contin¬ 
uum spectra with more than two independent frequen¬ 
cies between 178 and 8400 MHz over the whole sky. With 
these spectra it is now possible to assess the position and 
flux accuracy of the 19 large survey catalogues that we 
used for the spectrum determination (Table . Since the 
1400 MHz NVSS is a complete, deep, interferometric sur¬ 
vey, with by far the largest number (~ 1.8 million) of 
sources among these catalogues, we decided to use it as 
the reference for the positions of the other catalogues. 

2. Results 

To perform the cross-identifications (Vollmer et al. 2004) 
we used a least absolute deviation fit for the fit of the 
radio continuum spectra, which is more robust against 
outliers than a classical fit. This method has the dis¬ 
advantage, however, that the distribution of deviations is 
not Gaussian, but has in general a strong central peak 
and more extended wings. To circumvent this problem we 
made x^ fits instead of least absolute deviation fits to the 
frequency points found by SPECFIND, which results in a 
Gaussian distribution of deviations. 

The calculation of the offsets in the right ascension and 
declination of source positions is straightforward, and the 
uncertainty in the flux density is measured as follows: 


r. I 

‘-'extr 

where 5'cat is the flux density in the catalogue and S'extr 
is the flux density at the frequency of the catalogue as 
calculated using the fitted SPECFIND spectrum. Fig. ^ 
shows an example of the distribution of the differences 
in flux density (upper panel), of the offsets in right ascen¬ 
sion (middle panel) and of the offsets in declination (lower 
panel). The corresponding plots for the other catalogues 
can be found in the online version of this article. All three 
distributions can be fitted by a Gaussian of form 

iV = iVoexp( ^^~f°^' ) . (2) 

Since the accuracy of its position and the flux density 
of a radio source depend on the ratio of its flux density 
and the rms noise level of the survey, we calculate the 
uncertainties in the positions and flux accuracies for (i) all 
sources in a given catalogue, sources whose flux densities 
are (ii) smaller than 3 times the minimum flux density of 
sources in a given catalogue S'min, (iii) between 3 and 10 


Table 1. SPECFIND radio continuum source catalogue 
entries. 


name 

frequency resolution 
(MHz) (arcmin) 

Smin 

(mJy) 

number of 

sources 

Reference 

JVAS 

8400 

0.0055 

30 

2246 

(1) 

GB6 

4850 

3.5 

18 

75162 

(2) 

87GB 

4850 

3.5 

25 

54579 

(3) 

BWE 

4850 

3.5 

25 

53522 

(4) 

PMN 

4850 

3.5 

20 

50814 

(5) 

MITG 

4850 

2.8 

40 

24180 

(6) 

PKS 

2700 

8.0 

50 

8264 

(7) 

F3R 

2700 

4.3 

40 

6495 

(8) 

FIRST 

1400 

0.08333 

1 

811117 

(9) 

NVSS 

1400 

0.75 

2 

1773484 

(10) 

WB 

1400 

10. 

100 

31524 

(11) 

SUMSS 

843 

0.75 

8 

134870 

(12) 

B2 

408 

8.0 

250 

9929 

(13) 

B3 

408 

5.0 

100 

13340 

(14) 

MRC 

408 

3.0 

700 

12141 

(15) 

TXS 

365 

0.1 

250 

66841 

(16) 

WISH 

325 

0.9 

10 

90357 

(17) 

WENSS 

325 

0.9 

18 

229420 

(18) 

MIYUN 

232 

3.8 

100 

34426 

(19) 

4C 

178 

11.5 

2000 

4844 

(20) 


(1) Patnaik et al. (1992); Browne et al. (1998); 

Wilkinson et al. (1998) 

(2) Gregory et al. (1996) 

(3) Gregory & Gondon (1991) 

(4) Becker et al. (1991) 

(5) Wright et al. (1994; 1996); Griffith et al. (1994; 1995) 

(6) Bennett et al. (1986); Langston et al. (1990); 

Griffith et al. (1990; 1991) 

(7) Otrupcek & Wright (1991) 

(8) Fiirst et al. (1990) 

(9) White et al. (1998) 

(10) Condon et al. (1998) 

(11) White & Becker (1992) 

(12) Mauch et al. (2003) 

(13) Colla et al. (1970; 1972; 1973); Fanti et al. (1974) 

(14) Ficarra et al. (1985) 

(15) Large et al. (1991) 

(16) Donglas et al. (1996) 

(17) de Breuck et al. (2002) 

(18) Rengelink et al. (1997) 

(19) Zhang et al. (1997) 

(20) Pilkington & Scott (1965); Gower et al. (1967) 


times S'min and (iv) greater than 10 times S'min- For sources 
near the catalogue sensitivity the flux density uncertainty 
is typically ~ 20%, whereas at high signal to noise ratios 
the calibration uncertainty dominates. 

For source positions and flux densities the derived val¬ 
ues for Ao (the centre value, i.e., the mean systematic de¬ 
viation with respect to the reference, which is the NVSS 
for position and the value predicted by the composite 
spectrum from SPEGFIND for the flux density) and its 
standard deviation a are presented in Table |21 for all cata- 
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Fig. 1. Positional and flux density accuracy of the GB6 
survey with respect to the reference, which is the NVSS for 
position and the value predicted by the composite spec¬ 
trum from SPECFIND for the flux density. Upper panel: 
flux density (see Eq. middle panel: position offset in 
right ascension, lower panel: position offset in declination. 
Right hand panels: solid line: observations, dashed line: 
fitted Gaussian. The corresponding plots for the other cat¬ 
alogues can be found in the online version of this article. 


logues that were included in our analysis with SPECFIND 
(Vollmer et al. 2004), together with the accuracy in posi¬ 
tions and flux densities given by the catalogues’ authors. 
The columns are the following: (1) catalogue name, (2) Aq 
for the right ascension, (3) a for the right ascension, (4) 
Aq for the declination, (5) a for the declination (6) posi¬ 
tional accuracy according to the authors of the catalogue, 
(7) Aq for the flux density, (8) a for the flux density and 
(9) accuracy of the flux density according to the authors. 

3. Discussion 

We assumed a single or double power law for the radio 
spectrum in the spectrum-finding routine. Astrophysical 
sources show a single power law due to synchrotron emis¬ 
sion in a frequency range which is limited by the domina¬ 
tion of thermal emission by hot electrons at high frequen¬ 
cies (giving rise to a flat spectrum) and by synchrotron 
self absorption at low frequencies (giving rise to a flat¬ 
tening or inversion of the spectral slope). We take this 
into account by allowing a break in the radio spectrum 


if the object appears in a sufficient number of catalogues 
of different frequencies. For identification of two different 
slopes more than two independent frequency points must 
be associated with each slope, which is not the case for 
the majority of the sources. Therefore in general an over- 
or underestimate of the flux densities is expected calcu¬ 
lated using a single-slope SPEGFIND spectrum of sources 
with frequencies close to the edges of the frequency inter¬ 
val in which a single-slope radio continuum spectrum is 
expected. 

For the strong sources in the NVSS, we derived a po¬ 
sitional error of about 1” from our comparison with the 
sources in the 8400 MHz JVAS catalogue (Browne et al. 
1998), which has by far the highest positional accuracy 
(milliarcsec) among those used in our study. Note that the 
authors of the NVSS catalogue give positional uncertain¬ 
ties of < 1" for strong sources with flux densities greater 
than 15 mJy and of 1" for sources with flux densities at the 
survey limit. All positional uncertainties given in Table [21 
include this positional uncertainty of the NVSS. 

The scatter of NVSS flux densities with respect to 
those expected based on the SPECFIND spectra is about 
10% for flux densities greater than SS'min- The centre 
of the distribution of the JVAS flux density deviations 
is shifted to negative values, due to the insensitivity of 
the JVAS observations (made in the VLA A configura¬ 
tion, which lacks short baselines) to sources larger than 
7", despite the possible flattening of the spectra towards 
8400 MHz. The relatively large scatter of the JVAS flux 
density deviations (~ 18 %) is mainly caused by the VLA 
array parameters. 

The systematic offsets in right ascension and declina¬ 
tion are generally smaller than ~ 5 % of the survey resolu¬ 
tion for all catalogues with resolutions smaller than that 
of the NVSS (45"). Moreover, all offsets are smaller than 
the catalogues’ accuracy given by their authors. 

Since the 87GB is based on a subsample used for the 
GB6, which is more sensitive, we expect related results for 
their precisions in position and flux density. Indeed, the 
positional accuracies are comparable, and the precision 
in flux density is higher for the more sensitive catalogue 
(GB6), as expected. The positional uncertainties given in 
the GB6 and 87GB catalogues are for sources with flux 
densities greater than 10 times the minimum flux density 
of the catalogue. For the B2 survey, the accuracy in dec¬ 
lination is underestimated by ^ 50% by the authors. For 
the MIYUN catalogue we also find a lower positional ac¬ 
curacy (^ 20" in right ascension and declination for all 
sources and ^ 10" for sources with flux densities greater 
than 10 times the catalogue flux limit) than the > 5" given 
by the authors. 

There are no large systematic offsets in the flux densi¬ 
ties per catalogue. The flux density accuracy does not de¬ 
pend on the flux density of the sources for most of the cat¬ 
alogues (NVSS, JVAS, GB6, 87GB, BWE, PMN, MITG, 
PKS, FIRST, MRG, TXS, WISH, WENSS, MIYUN). The 
negative offsets of the 4C and MIYUN catalogues are most 
probably due to the flattening or turnover of the spectrum 
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Table 2. Positional and flux density uncertainties of selected radio continuum catalogues with respect to the reference, 
which is the NVSS for position and the value predicted by the composite spectrum from SPECFIND for the flux density. 


catalogue 

flux 

RA 

RA 

DEG 

DEC 

position 

Flux 

Flux 

Flux 


range 

centre 

G 

centre 

G 

accuracy 

centre 

G 

accuracy 


(xS'min) 

(arcsec) 

(arcsec) 

(arsec) 

(arcsec) 

(arcsec) 

% 

% 

% 

NVSS 

all 

< 3 
3-10 

- 

- 

- 

- 

< 5 

-0.9 

10.6 

3 


_ 

_ 

_ 

_ 


~ 0 

~ 10 




- 

- 

- 

- 


-0.9 

10.5 


JVAS 

all 

< 3 
3-10 

-0.1 


-0.2 

0.5 

< 1 

-11.1 

17.5 

3 


-0.1 

H 

-0.2 

0.4 


-12.6 

15.3 




-0.1 

■1 

-0.2 

0.6 


-7.2 

21.8 


GB6 

all 

-4.2 

17.9 

-3.9 

18.6 

8 

-1.3 

9.7 

10 


< 3 

-3.3 

23.5 

-1.0 

23.2 


-2.1 

9.3 



3-10 

-3.9 

14.4 

-4.2 

15.3 


0.3 

10.7 




-1.6 

10.0 

-1.5 

11.1 


-0.2 

10.2 


87GB 

all 

-3.2 

19.9 

-3.9 

22.9 

10 

8.5 

14.5 

14 


< 3 

-3.5 

24.6 

-3.7 

26.6 


9.0 

15.2 



3-10 

-1.7 

16.5 

-3.1 

19.7 


8.5 

14.1 




-1.4 

12.5 

-2.8 

16.7 


5.5 

11.9 


BWE 

all 

-1.6 

19.5 

-2.7 

22.4 

20 

-4.7 

11.7 

10 


< 3 

-1.7 

23.7 

-2.7 

25.6 


-6.0 

12.2 



3-10 

-0.8 

15.1 

-2.0 

18.4 


-2.8 

10.1 




-0.8 

12.7 

-1.9 

15.8 


-2.1 

10.9 


PMN 

all 

-15.3 

30.2 

-8.8 

22.7 

< 40 

-0.5 

5.4 

5 


< 3 

-3.5 

44.2 

-2.0 

34.1 


-0.5 

4.7 



3-10 

-14.9 

32.1 

-2.1 

22.9 


-0.3 

5.1 




-4.6 

15.2 

-8.1 

15.7 


-1.2 

8.9 


MITG 

all 

-3.1 

42.9 

-6.9 

44.0 

28 

-5.4 

20.4 

10 


< 3 

-0.8 

52.7 

-5.1 

55.2 


-8.6 

19.4 



3-10 

-1.7 

34.3 

-6.5 

33.5 


-2.0 

18.1 




5.7 

23.6 

-8.2 

23.8 


3.0 

19.1 


PKS 

all 

-1.6 

4.3 


2.9 

12 

-0.7 

12.6 

> 3 


< 3 

-1.0 

2.4 

BE 

2.1 


-6.1 

15.4 



3-10 

-1.1 

9.4 

-1.4 

8.6 


-1.1 

12.6 




-1.9 

1.9 

-0.8 

1.5 


0.7 

11.4 


F3R 

all 

1.4 

24.4 

-12.6 

27.6 

20 

3.1 

16.6 

20 


< 3 

1.4 

26.8 

-12.4 

30.1 


2.2 

17.9 



3-10 

2.0 

21.7 

-13.3 

24.0 


4.0 

15.0 




0.7 

18.5 

-9.9 

21.6 


4.5 

12.1 


FIRST 

all 

< 3 

-0.5 

0.6 

-0.2 

0.3 

< 1 

-2.5 

13.1 

5 



-0.5 

0.6 

-0.2 

0.3 


-2.5 

13.0 


WB 

all 

-0.1 

2.2 

-1.4 

2.2 

30 

0.9 

19.5 

15 


< 3 

-1.2 

2.2 

-2.1 

2.3 


-0.9 

20.7 



3-10 

-1.3 

1.4 

-1.5 

1.8 


2.7 

17.2 




-0.2 

1.2 

-0.7 

1.3 


4.8 

12.3 


SUMSS 

all 

< 3 

-0.4 

0.5 

-0.4 

0.3 

0.5 

2.0 

11.2 

3 



-0.4 

0.5 

-0.4 

0.3 


2.0 

11.1 


B2 

all 


12.5 

-4.1 

43.6 

10/20 

0.3 

14.5 

6 


< 3 


13.9 

-4.1 

53.6 


0.6 

15.4 



3-10 

0.9 

9.9 

-3.4 

30.7 


-0.4 

12.5 




-0.4 

8.7 

-1.3 

22.3 


2.7 

11.8 


B3 

all 

-1.5 

9.0 

-0.6 

14.1 

< 10 

-1.0 

10.3 

2 


< 3 

-1.2 

10.9 

-0.7 

19.3 


-0.2 

11.0 



3-10 

-0.5 

7.2 

-0.6 

10.5 


-2.3 

9.3 




-1.4 

5.2 

-0.2 

7.0 


-1.6 

7.9 
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Table 2. Uncertainties of catalogues with spectral indices (continued). 


catalogue 

flux 

RA 

RA 

DEG 

DEG 

position 

Flux 

Flux 

Flux 


range 

centre 

cr 

centre 

cr 

accuracy 

centre 

cr 

accuracy 


(xSmin) 

(arcsec) 

(arcsec) 

(arsec) 

(arcsec) 

(arcsec) 

% 

% 

% 

MRG 

all 

-1.3 

5.2 

-1.5 

6.3 

5 

1.7 

10.4 

7 


< 3 

-0.6 

5.4 

-1.3 

6.9 


1.7 

10.5 



3-10 

0.1 

3.4 

-1.0 

4.9 


1.6 

9.6 



> 10 

-0.9 

3.3 

-0.3 

3.8 


3.9 

11.7 


TXS 

all 

-2.1 

2.6 

-1.5 

2.3 

2 

2.1 

12.8 

5 


< 3 

-0.7 

2.0 

-1.5 

2.3 


2.2 

13.4 



3-10 

-1.0 

1.7 

-1.3 

1.7 


1.5 

11.8 



> 10 

-1.8 

2.2 

-0.7 

1.4 


3.3 

10.9 


WISH 

all 

< 3 
3-10 

-2.1 

2.0 

-2.9 

3.9 

> 2 

-4.7 

10.1 

10 


-0.2 

0.2 

-0.5 

0.9 


-0.2 

9.8 



> 10 

-2.1 

2.0 

-2.8 

4.0 


-5.0 

10.0 


WENSS 

all 

-1.7 

2.5 

-1.8 

2.4 

> 2 

0.7 

12.7 

6 


< 3 

-0.1 

0.3 

-0.4 

0.3 


-0.9 

10.9 



3-10 

-1.6 

2.2 

-1.4 

2.1 


1.0 

13.0 



> 10 

-1.8 

2.6 

-1.6 

2.5 


0.7 

12.7 


MIYUN 

all 

-1.7 

20.1 

-3.6 

25.9 

> 5 

-8.2 

17.2 

5 


< 3 

-1.4 

30.4 

-3.6 

40.9 


-7.5 

16.7 



3-10 

-1.2 

22.1 

-3.6 

30.2 


-9.8 

17.1 



> 10 

-1.7 

11.4 

-2.0 

13.3 


-6.1 

16.6 


4G 

all 

-2.3 

16.4 

-22.3 

173.4 

30/180 

-12.4 

15.4 

15 


< 3 

-1.9 

17.0 

-25.1 

189.3 


-12.3 

16.1 



3-10 

-1.7 

12.4 

-13.3 

96.2 


-12.4 

12.9 



> 10 

- 

- 

- 

- 


- 

- 



around 200 MHz. In general, the flux density accuracy 
lies in the range between 10% and 20% for all sources, 
and around 10% for sources with flux densities greater 
than 10 times the minimum flux density of the catalogue. 
This seems to indicate that the effective calibration er¬ 
rors are of the order of 10 % in all catalogues. The highest 
cr value for all sources are found for the JVAS, MITG, 
WB, and MIYUN catalogues. Our derived flux density 
accuracies (cr) are consistent with those given by the cat¬ 
alogues’ authors for about half of the catalogues (GB6, 
87GB, BWE, PMN, F3R, WB, MRC, WISH, 4G), whereas 
for the other catalogues we derived flux density accuracies 
that are more than a factor of two larger than those given 
by the authors. We did not detect a major problem in the 
flux density scale of any of the catalogues included in this 
study. 

This study shows that SPEGFIND results can be used 
efficiently to assess the position and flux accuracy of a 
large radio continuum catalogue. We therefore encour¬ 
age those responsible for future, as yet unpublished radio 
continuum catalogues to make use of our code to opti¬ 
mise the position and flux density calibration of their sur¬ 
veys. The SPEGFIND master catalogue, the manual and 
the SPEGFIND source files are publicly available in the 
VizieR database at the GDS^. 

^ Type SPEGFIND in the VizieR front page: 
http: / / vizier.u-strasbg.fr / viz-bin /VizieR 
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